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A next-to-leading order calculation for single top production including spin-dependent observables requires 
efficient techniques for the calculation of the relevant loop amplitudes. We discuss the adaption of dimensional 
regularization, the spinor helicity method and of tensor integral reduction algorithms to these needs. 



1. Single top production 

The top quark, discovered 1995 at Fermilab, 
is special among the quarks due to its large 
mass: Its lifetime is shorter than the character- 
istic hadronization time scale and therefore top 
bound states do not have time to form. From a 
calculational point of view this allows an imme- 
diate application of perturbation theory to top 
physics. Up to now top quarks have only been 
produced in pairs through the strong interaction. 
With the upcoming Run II at the Tcvatron and 
later with LHC one expects to produce a top 
quark also through an electroweak Wtb-vevtex. 
The principal production mechanisms for single 
top production are flavour excitation, W-gluon 
fusion, s-channel production and associated W- 
production. Some representative Feynman dia- 
grams are shown in fig. Q. Single top production 
will be an essential input for a direct measure- 
ment of the CKM matrix element Vtb- In addi- 
tion, a discovery of non-standard charged-current 
top couplings might give a hint on new physics. 
Since the electroweak decay of the top quark pro- 
ceeds so rapidly before any hadronization effects 
can take place, the decay products of the top 
quark are correlated with the top quark spin. Of 
particular importance will be the semileptonic de- 
cay t — > blv, since a detector signal corresponding 
to hadronic top decays will suffer from large QCD 
backgrounds. Single top production, in which a 
top quark is produced through a left-handed in- 
teraction, offers therefore an opportunity for po- 
larisation studies. Due to the short lifetime of 
the top quark, production and decay of a top 
quark should be considered together and the rele- 




s-channel Associated W-production 

Figure 1. Representative Feynman diagrams for 
single top production 



vant processes at the parton level arc for example 
u+b — > b+l+v+d for flavour excitation or u+g — > 
b+l+v+d+b for W-gluon fusion. Here it is 
understood that we may replace the (it, <i)-quark 
pair by (d, u), (c, s) or (s, c). Of course most di- 
agrams contributing to the two processes above 
do not contain a top quark at all. We performed 
a full leading-order calculation [ and compared 
the results to the top narrow width approxima- 
tion, to which only diagrams with an intermedi- 
ate top quark contribute. The top narrow width 
approximation significantly reduces the number 
of diagrams. Our results for the Tevatron are 
shown in table The first column gives the re- 
sults from the full calculation, in the second col- 
umn we required in addition that the decay prod- 
ucts of the top reconstruct to within 20 GeV to 
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Tevatron 


crtot 


±20 GeV 


narrow width 


Wg 


15.0 ± 0.4 fb 


14.3 ± 0.3 fb 


14.5 ±0.1 fb 


Wb 


87 ± 1 fb 


85 ± 2 fb 


87 ± 1 fb 




46 ± 1 fb 


32.3 ± 0.3 fb 


29.0 ± 0.2 fb 



Table 1 

Numerical results for Tevatron at 2 TeV for W-gluon fusion (Wg), flavour excitation (Wb) and s-channel 
production (qq). 



the top quark mass and the third column gives 
the results for the narrow width approximation. 
For W-gluon fusion we required three jets, two 
of them 6-tagged, for flavour-excitation two jets 
with one 6-tag and for the s-channel process we 
required two 6-tagged jets. Jets were defined by 
the hadronic /^-algorithm [ |) . This exclusive sig- 
nal definition allows us to distinguish the various 
processes. From table [j] we see that the narrow 
width approximation describes the cross section 
very well for W-gluon fusion and flavour excita- 
tion. The approximation is less satisfactory for 
the s-channel process. Here non-resonant terms 
seem to give a more sizeable contribution. We 
obtained similar results for the LHC. 



2. QCD corrections 

A leading-order calculation gives a rough de- 
scription of the process under consideration, 
but to reduce ambiguities due to the choice of 
renormalization or factorisation scales a next-to- 
leading order (NLO) calculation is required. Fur- 
thermore, if jets are defined an NLO-calculation 
models more accurately the internal structure of 
a jet. QCD corrections to the s-channel process 
and flavour excitation have been considered in [ 
||, f|. The calculation of QCD correction to W- 
gluon fusion is in progress [ |B| . This calculation 
should give a more reliable prediction on the px- 
spectrum of the 6-jet. Since the top narrow width 
approximation works well for W-gluon fusion, it 
is sufficient to calculate QCD corrections in this 
narrow width approximation. This simplifies the 
task considerably. The calculation of the relevant 
loop amplitudes can in principle be done entirely 
with conventional methods: One approach could 



be to calculate 

2 Re A Born Al 00 p (1) 

and to use the Passarino-Veltman algorithm for 
the reduction of tensor integrals [ |6|. One uses 
further the 't Hooft-Veltman prescription [ fjj for 
75, or, more efficiently, the reformulation accord- 
ing to S. Larin [ ||. Observables depending on 
spins can be treated within the spin density ma- 
trix formalism. This solves the problem of calcu- 
lating the one-loop amplitudes in principal, but 
leads in practice to complicated expressions. A 
more efficient approach is to calculate helicity am- 
plitudes and to square them numerically. Here 
the complexity grows linearly with the number 
of diagrams as opposed to a quadratic increase 
within the conventional approach. The calcula- 
tion of spin-dependent quantities is trivial, since 
helicity amplitudes carry the complete spin infor- 
mation. Efficient algorithm for the reduction of 
tensor loop integrals rely on Fierz and Schouten 
identities and "four-dimensional" regularization 
schemes like dimensional reduction [ ^| or the 
FDH-scheme [ |l0| are therefore favoured. (Of 
course these schemes are all variations of dimen- 
sional regularization.) Since weak interactions 
do not conserve parity, 75 makes it appearance 
and one carefully has to avoid inconsistencies in- 
herent in some of the four-dimensional schemes. 
We first review the spinor helicity method, before 
constructing a regularization scheme adapted to 
75, and comment finally on reduction algorithms 
for tensor integrals. 

3. Spinor helicity method 

The spinor helicity method expresses the polar- 
isation vectors for external gluons of momentum 
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k in terms of two-component Weyl spinors |p±) 
and (p ± | as 



S = 



(«-|7„l*-> 



(2) 



V2(qk) ' p \/2[feg] 

We have used the customary short-hand notation: 

(ij) = (Pi ~ \Pj+), [ij] = (Pi + \Pj~)- Hcre Q 
is an arbitrary light-like "reference momentum". 
The dependence on the choice of q drops out in 
gauge- invariant amplitudes. In the narrow width 
approximation we treat the massive top quark as 
an external state. For a massive spinor we use [ 



(3) 



Here, p is the four- vector of the massive fermion 




with p 



and po > 0, and q is an arbitrary 



null vector with qo > 0. It is easy to check that for 
these spinors the Dirac equations, orthogonality 
and completeness relations hold. It will be ad- 
vantageous to choose for q the momentum of the 
charged lepton from the top decay. This choice 
simplifies the factorisation formula in the narrow 
width approximation. 

4. Variants of dimensional regularization 

We first give an overview of existing dimen- 
sional regularization schemes. Naive dimensional 
regularization uses 

W,J"} = W / , {7 M ,75> = 0. (4) 

In the 't Hooft-Veltman scheme [ Q a D- 
dimensional object is split between a four- 
dimensional part and a (— 2e)-dimensional part: 

ki D) =k^+k^\ 7^^ 4) +7r e) -(5) 

Here 75 is defined as a generic four-dimensional 
object: 



(6) 



75 = l jf Vp ° iflfhfl^- 

As a consequence, 75 anticommutes with the first 
four Dirac matrices, but commutes with the re- 
maining ones: 



{7^,75} 



0. 



Ifj. i (5 



0. 



(7) 



Dimensional reduction [ |9j] continues the mo- 
menta to D < 4 dimensions, but keeps spinors 
and vector fields in four dimensions. 



2<$ 



(8) 



The D-dimensional metric tensor g^y acts also 
as a projection operator: 

9i D)p afJ = g$. 0) 

Naive dimensional regularization is not consistent 
In that scheme one can derive an equation like 



(D - 4)Tr 7 m 7 I /7 p 7<t75 = 0. 



(10) 



At D = 4 this equation permits the usual non- 
zero trace of 75 with four other Dirac matrices. 
However, for D ^ 4 we conclude that the trace 
equals zero, and there is no smooth limit D — > 4 
which reproduces the non-zero trace at D — 4. 
A similar equation can be derived in dimensional 
reduction, and dimensional reduction is therefore 
also algebraically inconsistent. Here the inconsis- 
tency can be related to the projection property 
eq. [)]. The 't Hooft-Veltman scheme is consistent. 
It should be mentioned that the 't Hooft-Veltman 
scheme violates the Ward identity for the axial 
current. This Ward identity relates the following 
three diagrams: 






The violation of the Ward identity is not a prob- 
lem per se, but one has to keep in mind, that one 
needs an additional finite renormalization to re- 
store the Ward identity. To prove the (algebraic) 
consistency of a regularization scheme there are 
two approaches: The axiomatic approach of Brc- 
itenlohner and Maison [ |l2| takes 7^, 75, g^ u etc. 
as a set of abstract symbols with a given set of 
relations and shows that these relations are self- 
consistent. The constructive approach of Wil- 
son and Collins [ |l3| tries to find a represena- 
tion for each object in some (finite or infinite- 
dimensional) vector space. Here algebraic con- 
sistency is ensured by the fact that one has an 
explicit representation. In the following we will 
follow the constructive approach [ p| . 
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5. Constructing a regularization scheme 

We start with a simple example. Suppose, we 
are given the natural numbers 1, 2, 3, ... together 
with the rule how to add two numbers and we 
are asked to construct or to define negative num- 
bers. We can do that as follows: We consider 
pairs of natural numbers and call two pairs, (p, q) 
and (k, I), equivalent, if 



p + l = k + q. 



(11) 



An addition of these pairs is defined by (p, q) + 
(k,l) = (p+k,q+l). The set of equivalence classes 
yields then the integer numbers. For example, the 
inverse of the class of (p, q) is the class of (q,p) 
and the neutral element is given by the class of 
(p,p). Given the set of integer numbers and a 
rule of multiplication we can repeat the exercise 
and construct the rational numbers. A further 
example would be to consider pairs (Vi , Vj ) and 
(Vk, Vi) from a set of vector spaces {Vi, Vi, V3, ...} 
and the equivalence relation 



Vi®V t ~ V k @V . 



(12) 



This construction is quite general and forms the 
basis of what mathematicians call K - theory. 
More formally, let A be an abelian semi-group. 
The if-functor associates to A an abelian group 
K(A), which is constructed as follows: Consider 
the equivalence relation on the set-theoretical 
product A x A. We put (a, b) ~ (a', b') when 
there exists a p € A such that 



a + b' +p = a' + b + p. 



(13) 



Then by definition K (A) = Ax A/ ~. We are now 
in the position to start the construction of our 
regularization scheme. Let V = {Vi, V2, Vi, ...} 
be a set of finite-dimensional vector spaces. V 
is an abelian semi-group with respect to the di- 
rect sum and the tensor product <g>. We use 
Grothendieck's K-functor twice to construct the 
corresponding abelian groups and define the rank 
of an quadruple [(Vi, Vj), (V k ,Vi)] by 



rank 



dim Vi — dim Vj 



(14) 



dim Vk — dim Vi ' 

Note the difference between the "rank" and the 
sum of dimensions dim Vi + dim Vj + dim 14 + 



dim Vi. We can construct a set of quadruples 
such that their ranks form a dense subset of K.. 
This construction can be extended to the com- 
plex case. Furthermore we can work with vector 
space of even dimension only. For a given quadru- 
ple of vector spaces of rank r and whose sum of 
dimension is 2m we define the integration by 

/ r/2-m f 
d r kf(k 2 ) = — / d 2m k 

OG 

■ Jdki(ki) r/2 - m -\f(e + ki) (15) 


The integration is well-defined, e.g. does not de- 
pend on the chosen representative. A different 
choice of a quadruple with the same rank r would 
yield the same result. In addition the definition 
satisfies the usual properties of an integration, 
e.g. linearity 

Jd r k{a.h + bf 2 ) - a Jd r k.h + b J d r kf 2 , (16) 
translation invariance 

J d r kf(k + q) = J d r kf(k), (17) 
the scaling law 

J d r kf(Xk) = \- r J d r kj(k) (18) 



and the normalization 

J cffcexp(-fc 2 ) = ir r/2 . 

To define the Dirac algebra we consider 



Va@V- 



2c, 



(19) 



(20) 



where V4 is a four-dimensional vector space and 
V-2e a quadruple of rank — 2e and total dimen- 
sion 2m. We recall that we can work with vector 
spaces of even dimension only. We now have two 
options to define the Dirac algebra. The first pos- 
sibility is to define Dirac matrices of dimensions 

2 2+m x 2 2+m Qver (y^ V-2e)- With 



75 



*7°7 1 7 2 7 3 



(21) 
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this yields the 't Hooft - Veltman scheme. The 
second option consists in defining Dirac matri- 
ces of dimensions 2 2 x 2 2 over V4 and of dimen- 
sions 2 m x 2 m over V~2e separately. This cor- 
responds to a four-dimensional scheme. By con- 
struction this scheme is free of algebraic incon- 
sistencies. Basically the distinction between the 
rank and the sum of dimensions avoids the pitfall 
of dimensional reduction. The scheme is speci- 
fied by a mapping of the Feynman rules from V4 
to V4© V-2e- The Dirac algebra in the numerator 
can effectively be performed in four dimensions. 
For the loop momentum appearing in the numer- 
ator of loop integrals one has 



- h 2 

— K (i-2e) 



h 2 



(22) 



The first term on the r.h.s can cancel a propa- 
gator, whereas the second term shifts effectively 
the rank of the loop integral from 4 — 2e to 

6 — 2e. As in the 't Hooft- Veltman scheme, the 
four-dimensional scheme violates Ward identities, 
which have to be restored by finite renormaliza- 
tions. 

6. Tensor integral reduction 

Loop integrals may involve the loop momen- 
tum in the numerator. To reduce these tensor in- 
tegrals to standard scalar integrals, several algo- 
rithm exists like the Pasarino- Veltman algorithm 
based on Lorentz invariance JJlq] , the Feynman 
parameter space technique [ |16[ or the partial in- 
tegration technique [ [l7|. These algorithm are 
general but have also some drawbacks: The first 
two introduce Gram determinants in the denomi- 
nator, whereas the last one might give rise to spu- 
rious poles in 1/e. For most practical applications 
more efficient algorithm exists. For the calcula- 
tion of helicity amplitudes loop momenta are usu- 
ally sandwiched between spinors. Consider the 
example 

W f h^Z. (23) 

7 7 J (2n)° (k 2 -m 2 1 )...((k-q n ) 2 -m 2 n ) { *> 

A tensor reduction according to Passarino- 
Vcltman would result in long intermediate expres- 
sions, but rewriting 



W*Mk = {k 2 -ml)+m\ 



(24) 



reduces the integral immediately. More system- 
atically, one uses the Fierz identity 

{A + \^\B+){C - \r\D-) = 2 [AD] (CB), (25) 

the Schouten identity for spinors 

{AS) {CD) = {AD){CB) + {AC){BD), (26) 

as well as the Schouten identity for four-vectors 



£(91,92,93,94)^ 
- 1 
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Q=l 



(k-q a )v° 



r 



abed 



s(fj,,q b ,q c ,q d ) 



(27) 



The algorithm proceeds in two steps: First one 
reduces tensor integrals to integrals with at most 
one power of the loop momentum in the numer- 
ator [ |l8| : If there at least two massless external 
legs, say p\ and p 2 , one writes 



1 



2pi ■ P2 



[(2fc-pa)j}i 



{2k ■ p\)i>2 



(28) 



The first two terms will cancel propagators, 
whereas a product of the last two terms can be 
rewritten as 

(l-|^|2-)<2-|^|l-) = (2*. Pl ) (2fe.ps) 
- {2 Pl ■ p 2 ) k 2 {i _ 2e) + (2pi • pi) kf_ 2e) , 

1 



(l-|fl2-)<2+|fll+> = 



<2- 



■[(2 Pl -p 3 )(2p 2 -k){l-\l 

- {2 Pl ■ k) {I - 
+ (2p 3 -fc)(l- 

-<l-|jMajM2- 



|3|1- 
12-) 



h 2 h 2 

ft (4-2e) K (-2e, 



(29) 



In a second step the remaining vector integrals are 
reduced [ |l9j. For pentagon or higher point in- 
tegrals one forms first the traces Tr±^i^ 2 ^3^4^5 
(here ± denotes the insertion of a helicity projec- 
tion operator) and uses then 

Tr± (k'p'ifofchfc) = -Tr T (tfifofofa) k 2 
-- (S ± e(k -Pi,P2,P3,Pa)) k 2 

+- (B 1 ± e(k,pi +p 2 ,P3,Pa)) k\ 
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- (B 2 ± e(k,pi,p 2 +P3,Pi}) k\ 



+ 2 ( B 3 ± £ ( fc 'Pl>P2,P3 +Pa)) k 3 
-- (B 4 ± e(k,pi,p2,pa)) k\ 

(30) 

The pi are the external momenta, the ki the mo- 
menta flowing through the loop propagators. For 
readability we dropped a subscript (4 — 2e) for 
the kf. Bi depends only on the pinched integral 
under consideration and is given for box integrals 
by 

Bi — st — m\m\ — m\rn\, (31) 

with s and t being Mandelstam variables and the 
nii denote the external masses. The derivation 
of eq. [50] is based on a combination of the spinor 
algebra approach with the dual vector approach 
[ [2(J . For box integrals one can use 
1 
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7. Summary 

Single top production is not only of interest for 
a direct measurement of Vtb or for a search for 
signals of new physics, but offers also the oppor- 
tunity of polarisation studies. In this article we 
discussed some technical issues relevant to a cor- 
responding NLO-calculation. 
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